Generating EPR beams in a cavity optomechanical system 



o 
o 



Oh; 

^3 ; 



> 
(N 
O 



oo 
O 



X 



Zhang-qi Yin^'^ and Y.-J Han^ 

^Department of Applied Physics, Xi'an Jiaotong University, Xi'an 710049, China 
'^FOCUS center and MCTP, Department of Physics, 
University of Michigan, Ann Arbor, Michigan 48109, USA 

We propose a scheme to produce continuous variable entanglement between phase-quadrature 
amplitudes of two light modes in an optomechanical system. For proper driving power and detun- 
ing, the entanglement is insensitive with bath temperature and Q of mechanical oscillator. Under 
realistic experimental conditions, we find that the entanglement could be very large even at room 
temperature. 
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Entanglement is the key resource of the field of quan- 
tum information. Light is the perfect medium to dis- 
tribute entanglement among distant parties. Entangled 
light with continuous variable (CV) entanglement be- 
tween phase-quadrature amplitudes of two light modes is 
widely used in teleportation, entanglement swap, dense 
coding, etc. [l| . This type of entangled state is also called 
Einstein-Podolsky-Rosen (EPR) state. The EPR beams 
have been generated experimentally by a nondegenerate 
optical parameter amplifier [2], or Kerr nonlinearity in 
an optical fiber Q. The later one is simpler and more 
reliable. The Kerr nonlinearity is used to generate two in- 
dependent squeezed beams. With interference at a beam 
splitter, the EPR entanglement is obtained between out- 
put beams. However, Kerr nonlinearity in fiber is very 
weak, which limits entanglement between output beams. 

It was found that strong Kerr nonlinearity appeared 
in an optomechanical system consisting of a cavity with 
a movable boundary [j, H, 0]. Besides, the single- mode 
squeezing could be made insensitive with thermal noise 
[5|, which makes the scheme very attractive. However, 
the frequency of output squeezed beams cannot be made 
identical, which makes interference difficult. Then it 
was generalized to two-mode schemes in order to gen- 
erate EPR beams without interference 0, H, [1, tlH ■ 
However, they are either very sensitive to thermal noises 
0) H, US) EH , or requiring ultrahigh mechanical oscillator 
Q ^ 10* to suppress thermal noise effects 0, which is 
2 to 3 orders higher than the present available parame- 
ters [13 • The practical scheme to generate EPR beams 
in an optomechanical system needs to overcome these 
problems. 

In this paper we propose a practical scheme to produce 
EPR beams in an optomechanical system, which consists 
of a whispering-gallery mode(WGM) cavity with a mov- 
able boundary. We find that, similarly as the single- mode 
scheme Q , the thermal noise in the two-mode scheme can 
be greatly suppressed by adiabatically eliminating an os- 
cillator mode. By precisely tuning the laser power and 
detuning, the oscillation mode is adiabatically eliminated 
and two output sideband modes are entangled. Unlike 
the cavity- free scheme Q , our scheme requires modest os- 
cillator Q. Besides, the output light is continuous in our 



scheme, other than pulse in Ref. 0. The most attractive 
feature of our scheme is that the entanglement between 
output beams is nearly not changed under different bath 
temperature and Q of the mechanical oscillator. Within 
the experimentally available parameters [isl . [l^ , we find 
the maximum two-mode squeezing could be higher than 
16 dB under room temperature. The entanglement of 
formation (EOF) between two modes is larger than 5 
(l5| . Since the coupling efficiency between cavity and 
fiber could be larger than 99% in the WGM cavity sys- 
tem [l^, we neglect the coupling induced noises in this 
paper. 




FIG. 1: (Color online) Experimental setup. Cavity 
modes a and b, which are driven by four lasers, couple 
to the mechanical mode am- 



As shown in Fig. [U we consider an optomechanical sys- 
tem consisting of a WGM cavity with a movable bound- 
ary. There are two cavity modes a and b with the same 
frequency but the opposite momentum. They are cou- 
pling with the same mechanical oscillation mode am and 
driven by four lasers, two from the right-hand side with 
frequencies ujl and ujl', the other two from the left-hand 
side with frequencies ujl and lul'. a-m and bin are input 
lights. a°^^ and a™* are output lights. Two lower mirrors 
have very high probability (> 99%) to reflect the driving 
lasers. So we neglect the reflecting induced noise for Uout 
and 6out ■ The system Hamiltonian is H = Hq + Hd -I- Hj , 
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where 



Hi —hv{a)b + ab^) + hrjiVmia^a + &^6)(am + o,]^) 



H.c. 



(1) 

Here a, 6 and the annihilation operators for the 

optical and mechanical modes, ujp and LUm are their angu- 
lar frequency. Qj with j — a,b is the driving amplitude 
and defined as flj = 2^ Pj /hujLT, where Pj is the input 
laser power and t = I/7 is the photon loss rate into the 
output modes, v is the coupling strength between cavity 
modes a and b. For the WGM cavity system, it ranges 
from 100 MHz to 10 GHz [H, H^. The dimensionless 
parameter 77 — {ujp/ujrn){xml R) is used to characterize 
optomechanical coupling, with Xm = y'h/mijjm the zero- 
point motion of the mechanical resonator mode [2l| . m 
its effective mass, and R a cavity radius. In typical WGM 
cavity systems we find 77 ~ lO"**. 

We define the normal modes oi = {a + b)/\/2 and 02 — 
{a-b)/V2. We suppose the conditions that Qa ~ = 
and rij, -f ri[, = are satisfied. The Hamiltonian can be 
written as 



H ^h{ujp + v)a\ai + h{ujp — i')a\a2 + hiOmal^am 



+ n(—aie ^ +—026 ^ 
+ hrjUJrnialai + a2a2)(al^ + am] 



H.c.) 



(2) 



where fii — + and = — We define the 
detuning Ai = lu]^ — Up — ly and A2 — uj']^ — ujp + v. As 
shown in Fig. [1] with beam splitters and Faraday rotator, 
we can get the output mode of a\ and a^- We assume 
both cavity and oscillator modes are weakly dissipating 
at rates 7 and 7™, respectively, where 7™ <C iOm- We can 
get quantum Langevin equations (2^ 

A / t \ ^7 7 / \ 



+V7< fori = 1,2, 



(3) 



a]aj ~ [iujm + — )am + V7^a^"(4) 

where thermal noise inputs are defined as corre- 
lation functions {a^^ {t),am{i')) — nmS{t — t'), 
{aZHt),a^Ht')) = {a^{t),a^{t')) = 0, 
{a^^{r),o}^{t')) = (aft(i),a;t(i,)) = (af (t), <(i')> = 
0, with n„i the thermal occupancy number of thermal 
bath for oscillator mode. We suppose cavity modes 
couple with vacuum bath. 

To simplify Eqs. ^ and we apply a shift to normal 
coordinate, aj — > aj + Uj, am am + P- otj and (3 are c 
numbers, which are chosen to cancel all c number terms 
in the transformed equations. We find they should fulfill 
the following requirements: (3 ~ — 77(|aip -I- |Q!2p), and 



i^ja-j+2iTfLOmaj{\ai\^-\-\a2\ 



2"J 



0. Because 



7m ^ uimi the imaginary part of (3 can be neglected. 
In the hmit Aj ^ 2ri'^Wm{\oti\'^ + |a2p), we find aj ~ 
Uj/ ^^'^ + 4A^. In the limit |a| ^ |(ap)|, the Langevin 
equations are linearized as 

hj = -iT]ujmajiam + a|„) -t- (zA' - -)aj -t- y/jaf[5) 
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(6) 

where j = 1,2 and A^- = A^ + 2T]'^ujm{\ai\'^ + |a2p)- We 
suppose A[ < and Aj > 0. We define S = (Aj - 
A'i)/2 - ujm and d = -{A[ + A^)/2. 

In the limit Wm S> 6, d^j,jm, the Langevin equations 
([5]) and ([6]) can be simplified as 



a 



0,1 — — idai 
02 = — ida2 
iSam - 



irjujmaiam - ^^i 



l-qUJmOi2K 
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V7<, 
V74", 

7»7 



/7mari 



(7) 

With proper detuning and input power, we can always 
tune the cavity mode amplitude ai = Q!2 = a. De- 
fine the Fourier components of the intracavity field by 
a{t) = -^j^^e-"^^^-^«^a{uj)duj. In the limit 5 > 
ti;,7m, we can adiabatically eliminate the am mode. We 
get a„i(ijj) ~ 77^(a*ai + aa\) — ^^^po'". Then we have 
quantum Langevin equations for ai(aj) and a\{—Lij) 

7 

—iujai{ijj) — — ig'ai{Lj) — igal2{—uj) cli{^) 



-iu}a\{—Lo) —ig'a\{—bj) + igai{uj) — ■^a2(— 
-.tin 



(8) 



where g = f]'^ \a\^u'^/5, 7^ = (ry|Q;|w™/(5)^7„i, and g' = 
,g + d. In Eq. ([8]), we neglect the phase of a because it is 
not important. 



Denote aiuj) 



and 



tion 
where 



We get the following matrix equa- 



Aa(cj) = V7«'"(w) 



/'""ml 



(9) 



A = 



«5 



Using boundary conditions a™* (w) = —a™{LLj) + ^aj{uj) 
for = 1,2, we can calculate output field as 



OUtf / 



[u) = G(c^)*4"^(c.) - iJ(c^)*a-(-c.) - /(c.)*a-(-c.), 

(10) 
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where G{u;) = (a;2 ^ ^ + - .g'^ - ig'j)/A{Lu), H{lu) = 
ig^/A{uj), I{uj) = {-iuj + ^ - ig' + ig)y/j^/ A{lj), and 

Let us define the dimensionless position and momen- 
tum operators of fields = [a1''\Lj) + a7*t(-w)] 

and P°"*(w) = [a™*(w) - a°"'^(-tj)]/i, for j = 1,2. We 
define the correlation matrix of the output field as Vij — 
((6?, +0C.)/2), where ^ = (Xr\ PP*, ^2°"*, ^2"')- We 
calculate the correlation matrix with Eq. (|10p . Up to 
local unitary transformation, the standard form of it is 



Vs 



fn fco; 

71 —kx 

kx n 

\0 -kx n 



(11) 



where n = {(w^ + 'J^ + ^ g'^f + [g'^ + ^2)^2 _,_ [(^ _^ 

.9' - 5)' + ^ ]77m(2n„ + l)}/|A(u;)|2, fc, = /^^T + l^, 
where = -2g-i{w^ + 1+9^- 9'^)l\A{u:)\\ ^24 = 

{2g'57^ + + 5' - 5)' + ^]77™(2n™ + 1)}/| A(c.)|2. 
This is the symmetric Gaussian state. The EOF for the 
symmetric Gaussian states is defined as [l5l | 



Ep = C+{n-kx) log2[C+(n - k^)] 
-C-{n - kx)log2[C-{n - k^)] 



(12) 



where C±{x) = {x"^^'^ ± x^^^Y /i. V describes an en- 
tanglement state if and only if n — k^ < 1. Based 
on the standard form of matrix (|lip . wc also find that 
+ X2)) = (-52(Pi - P2)) =n-kx. We define the 
two- mode squeezing as — — 101og]^Q(n — k^)- 




FIG. 2: (Color online) EOF for different temperature. 
\a\ is 1000, (5/27r = 10 MHz, d = O.O77. The maxima 
squeezing is larger than 16 dB when T = 300 K. 

We now estimate the bath noise influence in experi- 
mentally accessible conditions [l^. The cavity resonant 
frequency Wp = 27r x 300 THz. The oscillator frequency 
ojm = 2tt X 73.5 MHz. The mechanical Q factor is about 
30 000. The cavity and oscillator modes decay rates are 
7 = 27r X 3.2 MHz and 7™ = uJm/Q, respectively. The 



cavity radius is i? = 38 fim. The dimensionless coupling 
parameter is 77 ~ 10""*. We find if < c? ^ 7, the thermal 
noise does not decrease the maximum entanglement with 
strong enough input power. This is because we adiabat- 
ically eliminate the mechanical mode and suppress the 
effects of thermal noise. As show in Fig. [21 the tempera- 
ture change does not change the maximum entanglement 
with proper driving and detuning. But the higher the 
temperature, the less the entanglement spectrum width. 
The embeded figure shows that the maximum squeezing 
could be larger than 16 dB at room temperature. 




FIG. 3: (Color online) EOF for different cavity mode 
amplitude a. Here we adopt 5/27r = 10 MHz, T = 300 
K, 7p/27r = 3.2MHz, and d = O.O77. 



As shown in Fig. [31 the bigger the cavity mode am- 
plitude a, the larger the output entanglement. Because 

aj ~ rij I -I- 4 , the output entanglement is pro- 
portional to driving amplitude. But the peak of en- 
tanglement is splitted into two symmetric peaks when 
driving is very strong. The splitting distance is pro- 
portional to driving power. Increasing driving power 
can decrease the entanglement too. This is because 
adiabatical elimination condition w ^ (5 are not valid 
around peaks for very strong driving. So the driving 
power should be neither too big nor too small. For 
the specific a and 5, we find there is an optimum d 
which makes entanglement maximum and the entangle- 
ment peaks appear near w = 0. The optimum d is 
do = \/ {rf'Lo'^a^ / + 7^/4 — (r\uiOL)^ j i ^ corresponding 
to squeezing So = — 10 logjQ(4d2/72) and entanglement 
which is obtained from Eq. (|12p with n — kx = 4(^0/7)2. 
It is obvious that the higher the input power, the smaller 
the optimum d. In the mean time, we find that decreas- 
ing the mechanical Q factor nearly does not change the 
entanglement spectrum if d is around its optimum value 
and the condition ujm/Q ^ (5 is fulfilled. Leaving other 
parameters unchanged, Q could be as low as 300. Consid- 
ering the difficulty of increasing the mechanical oscillator 
Q, the above finding makes our scheme more practical. 
We also test the stability of our scheme. As shown 
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in Fig. [31 the optimum d is around O.O77 if a = 1000, 
6/2tt = 10 MHz. To maintain such high entanglement, 
we need to precisely control the d down to O.O27 ^ 27rx 60 
kHz. d is defined as d = -{A[ + A'^)/2 = -(Ai + A2)/2- 
Arj'^uj„i\ct\'^ ■ The higher entanglement is needed, the more 
precise detuning and driving power is required at the 
same time. To maintain the entanglement as high as Fig. 
m the laser spectrum width should be less than 60 kHz 
and the driving power fluctuation should be less than 1%. 
The lower entanglement between two beams is needed to 
maintain, the larger the optimum d is. Therefore higher 
fluctuations of detunina' and driving DOwer are allowed. 




FIG. 4: (Color onhne) EOF for different d. Here we 
adopt uJm/^TT = 73.5 MHz, T = 300K, 7™ = w^/30 000, 
7/27r = 3.2 MHz, and |a| = 1000. 



Before conclusion, we briefly discuss the approxima- 
tions we used. Our scheme needs the steady states ex- 



isting, which requires (ajaj) <C During numerical 

calculation, (a^a) is in the order of 10^, which is much 
less than |ap ~ 10^. The other two approximations are 
rotating wave approximation uJm ^ S, d, 7, 7m and adia- 
batical elimination 6 ^ co, 7m, which can be fulflUed inde- 
pendently. For a ~ 10'^, the driving amplitude fl is in the 
order of 10^^ Hz, which is much lower than the distance 
between adjacent cavity modes Auj — c/(i?no) ~ 5 x 10^^ 
Hz, where c is the light speed in a vacuum, hq the refrac- 
tive index of silica. Therefore the approximation that one 
laser only drives one cavity mode is valid. Laser power is 
needed in the order of 10 mW, which is available in the 
laboratory. 

In conclusion, we have proposed a scheme to gener- 
ate EPR lights in an optomechanical system. Two side- 
band modes, which couple with the mechanical mode, 
are driven by lasers. After adiabatically eliminating the 
the mechanical mode, we find that the output sideband 
modes are highly entangled. The higher power of the 
driving laser, the larger entanglement of the output light. 
To maintain the entanglement, we need to precisely con- 
trol the driving power and laser frequency at the same 
time. With proper parameters, the entanglement is in- 
sensitive to the thermal noise and mechanical Q factor. 
We test the scheme by experimental available parame- 
ters. Though in this paper we fucus on WGM cavity sys- 
tems, our scheme can be realized in other optomechani- 
cal systems, as long as the mechanical mode frequency is 
much larger than the cavity decay rate. 
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